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ABSTRACT 

We have computed the luminosity function for a sample of 389 field galaxies 
from the Canadian Network for Observational Cosmology cluster redshift survey 
(CNOCl), over the redshift range z = 0.2 — 0.6. We find Schechter parameters 
M* — 51og/i = —20.8 ± 0.4 and a = —1.3 ± 0.2 in rest-frame Gunn r, and 
^Bab -Slog/i = -19.6±0.3 and a = -0.9±0.2 in rest-frame Bab- We have also 
split our sample at the color of a redshifted but nonevolving Sbc galaxy, and find 
distinctly different luminosity functions for red and blue galaxies. Red galaxies 
have a shallow slope a ~ —0.4 and dominate the bright end of the luminosity 
function, while blue galaxies have a steep a ~ —1.4 and prevail at the faint end. 
Comparisons of the CNOCl results to analogous intermediate-redshift luminosity 
functions from the Canada-France (CFRS) and Autofib redshift surveys show 
broad agreement among these independent samples, but there are also significant 
differences which will require larger samples to resolve. Also, in CNOCl the red 
galaxy luminosity density stays about the same over the range z = 0.2 — 0.6, while 
the blue galaxy luminosity density increases steadily with redshift. These results 
are consistent with the trend of the luminosity density vs. redshift relations seen 
in the CFRS, though the normalizations of the luminosity densities appear to 
differ for blue galaxies. Comparison to the local luminosity function from the Las 
Campanas redshift survey (LCRS) shows that the luminosity density at 2 ~ 0.1 is 
only about half that seen at 2 ~ 0.4. A change in the luminosity function shape. 
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particularly at the faint end, appears to be required to match the CNOCl and 
LCRS luminosity functions, if galaxy evolution is the sole cause of the differences 
seen. However, it should be noted that the specific details of the construction of 
different surveys may complicate the comparison of results and so may need to 
be considered carefully. 



Subject headings: cosmology: observations — galaxies: luminosity function, mass 
function — galaxies: evolution 



1. Introduction 

The luminosity function of galaxies is a simple but fundamentally important quantity in 
the study of galaxy populations and their evolution. The luminosity function is of particular 
relevance to the problem of the excess counts of faint blue galaxies (e.g., [Koo fc Kron 1992| ; 



[Koo 199^ ) and will provide important constraints toward the resolution of this question. In 
particular, accurate determinations of the luminosity function at both low and high redshifts 
are crucial. Large wide-angle redshift surveys are providing precise measurements of the 
luminosity function in the local z ~ universe (e.g. |Lin et al. 1996a| ; [Marzke et al. 199^ ; 



[Loveday et al. 19^ ), a necessary baseline on which to anchor models of galaxy evolution back 



to higher redshifts. Also, recent smaller but deeper redshift surveys have provided direct 
measurements of the luminosity function to redshifts z 1 (e.g. [Lilly et al. 1995| ; [Ellis et 



al. 1996| ; [Cowie et al. 199(^ ; [Glazebrook et al. 1995[ ) , and have revealed clear evidence for the 



evolution of the luminosity function with lookback time. Moreover, that evolution depends 
strongly on galaxy type, such that the luminosity density of blue, star-forming galaxies 
appears to have increased substantially by z ~ 0.5, whereas that of red, more quiescent 
galaxies appears to have changed relatively little. 

In this paper we present the luminosity function for a sample of field galaxies obtained as 
part of the Canadian Network for Observational Cosmology cluster redshift survey (CNOCl). 
Although the CNOCl survey was optimized for obtaining cluster galaxy redshifts, a con- 
current sample of field galaxies was an important and necessary component needed in order 
to accomplish the main survey goal of an accurate measurment of VL from cluster dynamics 
( Parlberg et al. 1996| ). The 389 galaxies in the CNOCl field sample considered in this paper 
span the redshift range z = 0.2 — 0.6, and our sample size is comparable to that of other 
surveys at these intermediate redshifts. Moreover, available color information allows us to 
subdivide our sample by galaxy type. The aim of this paper is to compute the luminos- 
ity function for the CNOCl field galaxies and to compare our results to those from other 
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intermediate-redshift surveys. We describe the CNOCl data sample in § ^ and detail our 
methods in § ^ Our luminosity function results are then presented and discussed in § |[ We 
summarize our conclusions in § ^. 



2. The CNOCl Survey Data 

The CNOCl cluster redshift survey contains about 2600 velocities of cluster and field 
galaxies, observed in the fields of 16 high X-ray luminosity clusters spanning the redshift 
range z ~ 0.2 — 0.6. Photometric and spectroscopic observations were obtained using the 
Multi-Object Spectrograph (MOS) at the Canada-France-Hawaii Telescope (CFHT) during 
24 nights in 1993 and 1994. A detailed description of the observational and data reduction 
techniques is given in Yee, Ellingson, & Carlberg (1996; hereafter YEC); here we briefly 
describe some relevant details. 

The CNOCl field sample used in this paper is defined to lie within the redshift limits 
0.2 < z < 0.6, and outside the individual CNOCl cluster redshift limits given in Table 1 
of Carlberg et al. (1996). In addition, further redshift limits are introduced by the use of 
four band-limiting filters, designed to optimize galaxy spectroscopy for clusters in different 
redshift ranges (YEC). For each filter, we set the low redshift limit where [Oil] A3727 would 
be 50 A from the filter's blue end, and we set the high redshift limit where G-band would 
be 150 A from the filter's red end. (Note our high redshift filter limits are more conservative 
than the corresponding limits given in YEC, which were based on detectability of Ca H+K. 
The addition of G-band here and the resulting more stringent limits should minimize the 
effect of any possible incompleteness at the high- 2; ends of the filters, which we are still in- 
vestigating.) Selection of objects for spectroscopy was based on Gunn r magnitudes ( [I'huan 



|fc Gunn 19761 ). In order to optimize the number of cluster redshifts, the fraction of galaxies 
spectroscopically observed was designed to decline with apparent magnitude. Thus a mag- 
nitude selection weight Wm (where 1/wm is the fraction of objects with redshifts at a given 
apparent magnitude; see YEC) is assigned to each galaxy to properly weight its contribution 
in any statistical analyses (see below). In this paper we limit the sample to apparent magni- 
tudes for which Wm ^ 5. Gunn g magnitudes are also derived for each galaxy, and we further 
use a color selection weight wj^g — r) ( [Ellingson et al. 1996| ) to remove any residual color- 
dependent spectroscopic sampling effects not accounted for by the main magnitude selection 
weight. However, our results are little changed by inclusion of Wc- We use the observed g — r 
colors and interpolation among model galaxy spectral energy distributions ( [Uoleman et al. 



19801) to derive rest-frame colors and appropriate /c-corrections for determination of absolute 



magnitudes M, in both the Gunn r and Bab (|Oke 1972|) systems. (For this purpose, we also 
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need the transformations tab = r — 0.21 and qab = 9 + 0.05 given by Fukugita et al. 1995.) 
Throughout this paper we adopt a deceleration parameter go = 0.5 and a Hubble constant 
Hq = 100 h km s~^ Mpc~^; we use h = 1 if not otherwise indicated. Finally, we limit the 
absolute magnitude range to —22.0 < Mr < —17.5 or —21.5 < Mb^^ < —17.0, outside of 
which we have few observed galaxies. The final field sample consists of 389 galaxies. Table |l| 
summarizes the numbers, redshift limits, and apparent magnitude limits of our data sample. 



3. Methods 

We compute the luminosity function using two related methods which are unbiased by 
density inhomogeneities in the galaxy distribution. These are the parametric maximum- 
likelihood method of Sandage, Tammann, & Yahil (1979; hereafter STY) and the non- 
parametric stepwise maximum-likelihood (SWML) method of Efstathiou, Ellis, Sz Peterson 
(1988; hereafter EEP). We briefly describe these methods below. Fuller descriptions may be 
found in the STY and EEP papers, as well as in Lin et al. (1996a). 

Consider a galaxy i observed at redshift Zi in a fiux-limited redshift survey. Let mmin,j 
and minax.j denote the apparent magnitude limits of the field in which galaxy i is located, 
and also impose absolute magnitude limits Mi < M < M2 on the sample. Let 0(M) be the 
differential galaxy luminosity function which we want to determine. Then the probability 
that galaxy i has absolute magnitude Mi is given by 

Pi=p{Mi\zi)oc(j)iM,) / (l){M)dM , (1) 

/ Jma.x[M^iniz^),Mi] 

where Mmin(-2j) and Mmax(-2j) denote the absolute magnitude limits, at Zi, corresponding to 
the given apparent magnitude limits. We can then form a likelihood function, £, for having a 
survey of galaxies, with respective absolute magnitudes Mj, by multiplying the individual 
probabilities pi together, obtaining 

( /.min[M„,ax{2i)>A^2] 1 

In £ = V Hn 0(Mi) - In / (j){M)dM \ Wi + constant . (2) 

j=l \ Jmax[Mmin(2;,),Mi] J 

Here the weight Wi is a modification needed for the CNOCl survey in order to account for 
the apparent magnitude and color selection effects described in § ^ Wi = Wm,iWc^i (cf. |Zucca 
8t al. 19941 ). 

In the STY method one assumes a parametric model for 0(M), and the parameters 
describing 0(M) are determined by maximizing the likelihood £, or equivalently ln£, with 
respect to those parameters. In our case we take as our model for 0(M) the Schechter 
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function ( |Schechter 1976| ) 

0(M) = (0.41nl0) 0* [ioO-4(M*-A/)]i+a exp[-10°-<^^'-*^)] , (3) 

and use the STY method to find the characteristic magnitude M* and the faint-end slope 
a. We emphasize that in the STY method, the determination of M* and a are unbiased 
by the presence of galaxy density fiuctuations in the survey, unlike the case for traditional 
least-squares luminosity function estimators (e.g., [Felten 1977|) . On the other hand, the 
normalization (p* drops out in equation (|l|) and has to be determined separately as described 
below. Error ellipses in the M*-a plane may be drawn by finding the contour corresponding 
to 

ln£ = ln£max-|Ax^ , (4) 

where Ax^ is the change in appropriate for the desired confidence level and for a 
distribution with 2 degrees of freedom (e.g., Ax^ = 6.17 at 2a confidence). 

Alternatively, one does not have to assume a particular functional form for 0(M). 
Rather, in the EEP stepwise maximum likelihood method (SWML), the luminosity func- 
tion is taken to be a series of Np steps, each of width AM in absolute magnitude: 

0(M) = 0fc , Mfc - AM/2 <M <Mk + AM/2 , k = 1, . . . , Np . (5) 

Here the likelihood is maximized with respect to the (pk, which are solved for via an iterative 
procedure. Also, we can estimate the variances of the (pk by calculating the diagonal elements 
of an appropriate covariance matrix. See EEP for details of the procedures used to estimate 
the 0fc and their errors. 

To calculate the normalization 0* of the Schechter function (^, as well as the mean 
galaxy number density p, we do the following. For galaxies i within redshift limits zi < Zi < 
Z2 and absolute magnitude limits Mi < Mj < M2, we estimate p by 

P= y , (6) 

where V is the appropriate (comoving) volume, and S{z) is the selection function defined by 
S{z) = / 4>{M)dM / / (P{M)dM . (7) 



/max[Afn,in(^),Mi] / JA/i 

The Schechter function normalization 0* is then just 

6* = P 

J^,[^ <P'{M)dM 



(8) 
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where 0' is the Schechter function with 0* set to one. We use bootstrap resamphng (e.g., 
Barrow et al. 198^) to estimate the errors on p and 0* contributed by uncertainties from 



fitting M* and a, and from samphng and weighting fluctuations. This method does not 
include errors arising from galaxy density fluctuations, which we instead estimate using an 
appropriate integral over the galaxy clustering power spectrum P{k): 



{6-p/py = j^^Jd'kP{k) \W{k)\\ (9) 

Here V denotes the volume of the particular sample under consideration, and W(k) is the 
window function for that volume, deflned by 

W{k) = ^ j^dh e'^'' . (10) 

(That is, in V we explicitly account for the "pencil-beam" geometry of the CNOCl flelds.) 
For P{k) we adopt the local result derived from a 19000-galaxy sample drawn from the 
Las Campanas Redshift Survey (LCRS; |Lin et al. 1996bD , but multiplied by 1/(1 + z)'^ to 
account for linear evolution (e.g., [Peebles 1993| , pp. 528-529) at the higher redshifts sampled 



in CNOCl. The flnal errors on 0* and p consist of the quadrature sum of the bootstrap 
resampling and density fluctuation errors; it turns out for our sample that these two sources 
of error make roughly equal contributions. We will also calculate luminosity densities pi for 
our galaxies, using 

EiW, io-o-4^^V^l(^. 



where 



PL = y , (11) 

Sl{z) = / 10-°-^^V(M)dM / / 10-°-^*V(M)dM ; (12) 

Jmax[Mmin(^),Mi] / J -oo 

that is, we sum over the luminosities of our observed galaxies, but weighted by the factor 
Sl{z), which uses the luminosity function to extrapolate for the luminosity of unobserved 
galaxies lying outside the accessible survey flux limits. Note we use a flnite Mfamt (—17-1- 
51og/i for both r and Bab bands) instead of +oo to prevent the denominator of Sl{z) from 
blowing up when a < —2, as sometimes occurs during bootstrap resamplings of our blue 
galaxy subsample (§ Also, note that our survey apparent flux limits prevent us from 
measuring the luminosity function fainter than our chosen Mfamt (cf. end of § |^). We 
estimate errors on pi with the same method used on 0* and p. 

We have thus far neglected the effects of photometric errors. This is actually a reasonable 
simpliflcation, as 99% of our sample galaxies have estimated errors of < 0.1 mag. Assuming 
a gaussian magnitude error distribution with dispersion a = 0.1 mag, we can correct for 
the effects of photometric errors on the luminosity function using the method described in 
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EEP § 3.5. We find that neglecting photometric errors only biases M* and a by at most 
AM* = —0.02 and Aa = —0.01, much smaller than our la uncertainties (Table 0). We can 
thus safely neglect the effects of photometric errors in the rest of this paper. 

Finally, note that for galaxies within our adopted filter redshift limits, the magnitude 
selection weights Wm will be slight overestimates, for the following reason. Recall that Wm 
is basically the ratio (at a given apparent magnitude m) of the total number of galaxies to 
the number of galaxies with redshifts. However, since we include galaxies both inside and 
outside the filter redshift limits when computing this ratio, and since the redshift success rate 
will be lower for galaxies outside the filter limits than for those inside, Wm will be overesti- 
mated for those galaxies inside the limits. The amount of overestimate can be approximately 
calculated once some simplifying assumptions are made. We adopt specific values for M* 
and a (those we find below in § ^; we further assume that the luminosity function does not 
change with redshift), and we take our redshift sampling rate to be some constant within a 
given set of filter redshift limits. It turns out that the effect is approximately independent of 
apparent magnitude, so that the luminosity function shape is not affected. The luminosity 
function normalization will, however, be too large by roughly 20± 10% for our adopted filter 
limits, where the uncertainty represents the range of scatter observed for different apparent 
magnitudes and filters. We do not explicitly include this systematic correction in our lumi- 
nosity function normalizations, as it is comparable to our random errors, as well as somewhat 
complicated to deal with exactly (e.g., we need to model the change of luminosity function 
with redshift, as well as iterate the luminosity function and "renormalization" calculations 
until convergence). However, we will need to keep the approximate 20% correction in mind 
when comparing our results to those of other surveys, as we do below. 



4. Results 

Our luminosity function fit results are given in Table H, for both rest-frame Gunn r and 
rest-frame Bab- Figures |I| and ^ show the 2a error ellipses in the M* and a parameters, 
and Figures |^ and ^ plot the actual luminosity functions 0(M). The results are shown 
for the full 389 galaxy sample, as well as for two subsamples of galaxies whose rest-frame 
colors are either redder (209 galaxies) or bluer (180 galaxies) than that of an Sbc galaxy 
(model of [Coleman et al. 1980| ). Overall, the full sample luminosity functions have Schechter 
parameters M* = —20.8 and a = —1.3 in Gunn r, and M* = —19.6 and a = —0.9 in Bab- 
Note the clear distinction between the red- and blue-subsample luminosity functions shown 
in Figures § and |^: the red subsample has a shallow a = —0.4 and dominates the galaxy 
population at the bright end {Mr,MB^g < —19.5), while the blue subsample has a much 
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steeper a = —1.4 to —1.5, and consequently dominates at the faint end (M^, Ms^g > —19.5). 
The error eUipses for the red and blue subsamples are clearly separated for both Gunn r and 
Bab, as seen in Figures ^ and 0, respectively. 

We next compare our results against those from the two other largest surveys with 
5-band luminosity functions available at comparable intermediate redshifts, specifically the 
Canada-France Redshift Survey (CFRS; [Lilly et al. 1995| ) and the Autofib Redshift Survey 
( piis et al. 19961 ). We start with the CFRS. In Figures | and | we plot the CFRS ^As-band 
results for a sample of 208 galaxies in their z = 0.2 — 0.5 bin (their "best" estimates). The 
CFRS sample has also been divided into red and blue subsamples at the same Sbc color cut 
used in CNOCl, and we show the corresponding results as well. Note first from Figure ^ 
that the CFRS M* and a values are just outside the respective CNOCl 2a error ellipses. 
Remembering that the corresponding CFRS error ellipses will be larger (the CFRS sample 
is about half the size of CNOCl), we find that the CNOCl and CFRS M* and a values 
agree at better than the 2a level for each of the all, red, and blue samples. However, when 
we plot the actual luminosity functions in Figure ^, we find that although the shapes (i.e., 
a) of corresponding CNOCl and CFRS luminosity functions agree fairly well, the CNOCl 
results tend to show brighter M* values and higher normalizations, especially in the case 
of the blue subsample. Recall that the CNOCl normalization is likely an upper bound and 
may have to be decreased by about 20%, which would help improve the match with CFRS. 
An additional decrease in the CNOCl M* by 0.2 mag (cf. Figure ^ would be needed to 
bring the CNOCl and CFRS full-sample luminosity functions into very close agreement. 

We consider next some possible causes for the differences, though none of them are 
entirely satisfactory. Extinction, which has been neglected in our analysis, should only 
have a small effect (as our photometry was in the r-band) and would actually brighten our 
M* values. Large-scale structure may bias the CNOCl normalization high, as our field 
samples are in the neighborhood of rich clusters. However, we find that our results are little 
changed even when we triple the adopted cluster redshift limits (from about ±3000 km s~^ 
to ±9000 km s~^ from the cluster centers) in order to reduce the effects from potential 
large-scale structure associated with the clusters. Moreover, one expects the discrepancies 
to be larger for red galaxies (as they are more common in dense environments), but the 
observed differences are actually larger for blue galaxies. Also, the choice of spectral energy 
distributions used here to calculate fc-corrections and to make the red/blue cuts is the same 
as that in CFRS. We also verified that we could reproduce the CFRS luminosity function 
results using our code on their redshift catalog data. A final possibility is that there remains 
some unaccounted systematic offset of order 0.2 mag in the photometry zero-points, arising 
from the completely independent ways the photometry was obtained and reduced in the 
two surveys. In the end, the differences may simply reflect actual sampling fluctuations. It 
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is nevertheless encouraging that the luminosity function shapes of corresponding CNOCl 
and CFRS samples are quite similar. Also, even before the 20% correction to the CNOCl 
normalization, the luminosity densities from the two surveys do not differ at more than the 
2(7 level (see also below). Overall the CNOCl and CFRS results are roughly consistent, 
showing the same luminosity function shapes, but with normalization and M* differences, 
both at about the 20% level. We note that the CNOCl sample provides somewhat better 
luminosity function measurements in the range z = 0.2 — 0.6 because of its larger sample 
size. On the other hand, the full CFRS sample extends to 2; ~ 1 and provides a longer 
baseline for studying the evolution of red and blue galaxy luminosity functions, a point we 
return to below. 

Next, in Figures |^ and |^ we compare our results against the z = 0.15 — 0.35 and 
z = 0.35 — 0.75 luminosity functions from the Autofib survey. (The Autofib survey results 
are in the bj system, which is close enough to Bab that we will not apply any zero-point 



corrections; cf. [Fukugita et al. 1995| . Also, there are no results for red and blue subsamples.) 
Here the agreement is worse than it was with CFRS. Both Autofib samples show steeper faint- 
end slopes than those of CNOCl and CFRS. The M* and a values are discrepant at the 2 to 
3(T level compared to CNOCl (the corresponding error ellipse for the Autofib z = 0.15 — 0.35 
sample is smaller than that of the CNOCl all sample, while that for z = 0.35 — 0.75 is 



comparable in size to that of CNOCl; see Figure 11 of |Fllis et al. 1996 ). Examination of 



Figure ^ shows that the CNOCl, Autofib, and CFRS luminosity functions are only broadly 
consistent with each other, with Autofib samples showing the steepest slopes and lowest 
normalizations. It is curious that the Autofib results actually resemble that of our blue 
subsample well (though this is not to suggest that they are incomplete in red galaxies). 
Some of the same points raised above in discussing the CNOCl and CFRS differences apply 
here as well. Also, though the Autofib sample is the largest among the three surveys, 
there are some differences in their survey construction compared to CNOCl and CFRS. 
Unlike CNOCl or CFRS, the Autofib sample is actually a combination of a number of 
disparate redshift surveys, carried out with different instruments and telescopes. Instead of 
CCD photometry, the Autofib survey is based on photographic plate photometry. Also, the 
original photometry is in the bj band, rather than at longer wavelengths like r (CNOCl) or 
/ (CFRS). Consequently the /c-corrections needed to convert to rest-frame B are largest for 
the Autofib samples. Ellis et al. (1996) and references therein detail a careful treatment of 
these various issues. While we do not claim that any of these is the cause of the luminosity 
function differences, we do suggest that the details of the construction and selection of galaxy 
samples from one survey to another are complicated enough that they may systematically 
affect comparison of results among different samples. The less than ideal agreement shown 
in Figure ^ indicates that even larger samples, with well-defined selection criteria, are needed 
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to make a definitive determination of tlie luminosity function at intermediate redshifts. 

Now, we proceed to examine tlie issue of galaxy evolution within the CNOCl sample. 
We have divided the CNOCl samples by redshift to look for changes in M* and a as a 
function of z. We find there to be no significant changes in M* or a with redshift, within 
the respective 2a error ellipses of the 2;-divided samples. This holds for the full sample, as 
well as for the red and blue galaxy subsamples. However, because the error ellipses do get 
quite large when the redshift cuts are applied, we are not able to rule out (or in) changes in 
M* of order 0.5 magnitude over 0.2 < z < 0.6. 

We will instead examine changes in the galaxy luminosity density p^, as a function of 
2;, and adopt the simple procedure of fixing M* and a at the values given in Table || for 
all z. We plot the S^s-band luminosity density vs. z in four redshift bins in Figure 
for the all, red, and blue samples. (To facilitate comparison with other surveys, we have 
rescaled our pi, results for Mfamt = —17 + Slog/i to Mfaint = +C)o; recall equation [|12|-) The 
errors plotted are la estimates as described in § |^, with roughly equal contributions from 
bootstrap resampling and density fluctuation uncertainties. We see from the figure that p^ 
for the full sample rises at higher 2;, and that this increase is caused predominantly by the 
blue galaxy population. The red sample p^ does not exhibit a trend with redshift, while the 
blue sample pi, rises steadily with z. If we use just two bins (for simplicity), z = 0.2 — 0.4 and 
z = 0.4 — 0.6, the red galaxy p^ in the high- 2; bin is just 10% more than that in the low- 2; bin, 
while the blue galaxy nearly triples, albeit with large errors. These results may be directly 
compared to the luminosity density vs. redshift results derived from the CFRS ( [Lilly et"ar| 
|1996|) which we also show in Figure ^ We find that the shapes of the CNOCl pl vs. z relation 
are fairly similar to the CFRS results, which were derived over a longer redshift baseline. 
However, the CNOCl luminosity densities appear systematically higher, especially for the 
blue sample. If we apply the likely 20% normalization reduction for CNOCl discussed before, 
the results agree better, as also shown in Figure ^. In any event, the CNOCl trends are at 
least qualitatively consistent with those seen in the CFRS: strong evolution of the luminosity 
density of blue galaxies, but little evolution of that of red galaxies. Moreover, Schade et al. 
(1996) have found that the surface brightness of disk galaxies in the CNOCl survey increases 
with z in a way consistent with the shape of the CFRS pi-z relation; our current results 
thus also corroborate the earlier CNOCl findings. That the differences between CNOCl 
and CFRS lie in the absolute normalizations rather than in the shapes suggests that within 
each survey, we are basically measuring the same galaxy luminosity density trends, but that 
there is some unaccounted systematic normalization or scaling difference, particularly for 
blue galaxies, which causes the disagreement when we compare across the two surveys. 

We have also checked that our red/blue luminosity function and luminosity density 



evolution differences are robust to potential errors in assigning galaxies to the red and blue 
samples (e.g., a potential wavelength-dependent error in our model Sbc galaxy spectrum 
can lead to a redshift-dependent type assignment error). We do this check by more finely 
dividing our sample into three roughly equal parts by color. We find that the steep faint 
end of the blue luminosity function and the increase of the blue luminosity density at higher 
redshift both result primarily from galaxies in the bluest third of our sample, rather than 
from blue galaxies near the Sbc dividing line. In the latter case our results would have 
been more sensitive to potential type-assignment errors, but as the former case is the actual 
situation, our conclusions regarding the red/blue differences should not be affected. 

Finally, we compare our r-band results against those from the local sample of 18,678 
Las Campanas Redshift Survey (LCRS) galaxies, the largest galaxy sample for which the 
luminosity function has been computed ( |Lin et al. 1996a| ). The LCRS galaxies have an 
average redshift z = 0.1, and are also observed in the red0. From Figures |l| and |^, we find 
that the local LCRS luminosity function is clearly different from the intermediate redshift 
CNOCl results. The CNOCl results show a brighter M*, a steeper faint-end slope a, as 
well as a higher normalization compared to the LCRS. We note that if the differences result 
purely from galaxy evolution, that evolution needs to to be luminosity dependent as well, 
as the differences become larger at the faint end of the luminosity function. In terms of 
luminosity densities, the LCRS has pi = (1.9 ± 0.2) x 10* Lq h Mpc~^, about half that 
of the CNOCl sample. The increase in luminosity density at intermediate redshifts relative 
to the local value appears to be a robust conclusion, as the CNOCl results are corrobo- 
rated by independent samples hke CFRS and Autofib (despite the differences among the 
intermediate-redshift samples). However, a careful consideration and account of systematic 
effects will be especially important for a definitive assessment of the causes of the differences. 
For example, different surface brightness limits of local vs. deeper surveys ( [Ferguson fc Mc-] 
Gaugh 1995| ), or different photometry methods (different aperture sizes, limiting isophotes, 
etc.) can cause differences in the luminosity function that have nothing to do with actual 
galaxy evolution. Though beyond the scope of this paper, such a detailed modelling of both 
evolutionary and observational effects is planned for the ongoing CN0C2 redshift survey, 
a successor of CNOCl, which should eventually provide a much larger sample of several 
thousand intermediate redshift field galaxies for study. 



''Specifically an isophotal, "hybrid" Kron-Cousins i?-band (details in Tucker 1994 and Lin et al. 1996a), 
whose zero-point is within « 0.1 — 0.2 mag of the Gunn-r system used for CNOCl. 
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5. Conclusions 

We have computed the rest-frame Gunn r and Bab luminosity functions for a sample 
of 389 field galaxies from the CNOCl redshift survey, over the redshift range z = 0.2 — 0.6. 
We find Schechter parameters M* — Slog/i = —20.8 ± 0.4 and a = —1.3 ± 0.2 in Gunn 
r, and M^^^ — Slog/i = —19.6 ± 0.3 and a = —0.9 ± 0.2 in Bab- Samples of red and 
blue galaxies, cut at the rest-frame g — r color of an Sbc galaxy, show distinctly different 
luminosity functions. Red galaxies have a shallow slope a ~ —0.4 and dominate the bright 
end of the luminosity function, while blue galaxies have a steep a — 1.4 and prevail at the 
faint end. Table |^ summarizes our results. 

Comparisons of the CNOCl results to independently-determined intermediate-redshift 
luminosity functions show broad agreement with results from the CFRS and Autofib redshift 
surveys. However, there are questions about particular differences in the luminosity functions 
which will require larger samples to resolve. Calculation of the -Byi^-band luminosity density 
shows that the CNOCl red galaxy luminosity density is about the same from z = 0.2 — 0.6, 
but that the blue luminosity density over z = 0.4 — 0.6 is nearly three times that for z = 
0.2 — 0.4. These trends are consistent with those derived from the CFRS (see Figure |^), 
except for a normalization difference resulting primarily from blue galaxies, and are also 
consistent with the redshift evolution of disk galaxy surface brightness observed in CNOCl. 
Comparison to the local luminosity function from the Las Campanas redshift survey implies 
that the luminosity density at 2; ~ 0.1 is only about half that seen at 2 ~ 0.4. Also, 
luminosity- dependent evolution, increasing at the faint end, would seem to be required to 
match the CNOCl and LCRS luminosity functions, if galaxy evolution is the sole cause of 
the differences seen. 

However, an underlying caveat throughout our comparisons of CNOCl results with 
those of other surveys is that the particular details of the construction of different surveys 
may have an important effect on how we interpret the results. Along these lines, we have 
pointed out some potential systematic causes of the differences among the various survey 
luminosity functions we considered, though we have not attempted a detailed investigation. 
It is probably fair to say that control of systematic effects within a single survey is simpler 
than across different surveys, and perhaps that is why the red/blue luminosity density trends 
appear to have the same shapes but different normalizations in CNOCl and CFRS. In any 
event, the CNOCl field sample is still relatively small, and a larger sample will definitely 
help confirm or reject some of the trends that we have seen in CNOCl, as well as help resolve 
the differences seen vs. other surveys. In particular, the CN0C2 field redshift survey should 
provide us with an order of magnitude larger sample of several thousand objects to examine 
the luminosity and other properties of galaxies at intermediate redshifts. We thus anticipate 
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a similar but more detailed analysis for the CN0C2 data in the near future, and a thorough 
investigation of the evolution of galaxies from intermediate to low redshifts. 

We thank the many CNOC team members who participated in the observation and 
reduction of data for the survey. We thank Simon Lilly and David Schade for helpful discus- 
sions and for providing spectral energy distributions and related information. Wc also thank 
the referee David Koo for helpful suggestions. Financial support from NSERC and NRC of 
Canada is also gratefully acknowledged. 
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Table 1. Sample Parameters 



Filter^ Redshift Limits Apparent Magnitude Limits ^ 



Z2 99 

Z3 45 

Z4 217 

Z5 27 



0.20-0.27 
0.27-0.38 
0.21^- 0.43 
0.45 - 0.60 



18.0-20.5 
19.0 - 21.5 
19.0-21.7 
20.0 - 22.0 



^Names of band limiting filters used for spectroscopy; see text 
and YEC for further details. 

•^The samples are further restricted to the Gunn r absolute mag- 
nitude range —22.0 < Mr — 5 log h < —17.5. The numbers remain 
the same for the alternative restriction —21.5 < Mb^b < —17.0, 
except that there is one additional galaxy in the Z4 sample. 

'^Gunn r magnitudes. 

^The value given in YEC was 0.27 and was based only on de- 

tectability of [Oil] A3727. The present value is extended down 
to 0.21 based on the added dctectability of [OIII] AA5007,4959 
at those low redshifts where [Oil] A3727 would disappear off the 
blue end of a spectrum. Experience has shown that in the CNOCl 
sample the [Oil] and [OIII] features are always detected simulta- 
neously whenever they are all potentially visible. 



-16- 



B 



Table 2. Luminosity Function Parameters 



Sample N M*-51og/i^ a 0* '^'^ Pi ^'^ 
Gunn r 

All 388 -20.80 ±0.40 -1.25 ±0.19 0.021 ± 0.008 4.3 ± 0.5 

Red 209 -20.24 ±0.30 -0.42 ± 0.28 0.024 ±0.005 2.3 ±0.3 

Blue 179 -20.24 ±0.52 -1.47 ±0.32 0.013 ± 0.009 1.8 ±0.4 



AB 



All 389 -19.63 ±0.25 -0.89 ± 0.21 0.042 ±0.010 3.6 ± 0.4 
Red 209 -19.46 ± 0.30 -0.38 ± 0.29 0.024 ± 0.005 1.8 ± 0.3 
Blue 180 -19.85 ±0.50 -1.44 ±0.32 0.014 ± 0.009 2.0 ±0.4 



^Absolute magnitude restrictions —22.0 < — 51og/i < —17.5 for 
Gunn r, and —21.5 < Mb^b — 51og/i < —17.0 for Bab-, are applied in the 
definitions of the samples and the calculations of the luminosity function 
and associated parameters, except that pl is computed for — oo < M — 
51og/i < -17. 

''Errors are la one-parameter errors determined from the STY fits. See 
Figures 1 and 2 for the full ^u;o-parameter M*-a error eUipses. 

'^Units are Mpc~^ mag~^. Errors are la and are estimated as described 
in the text. 

"^Units are 10^ Lq h Mpc"^. We take Gunn tq = 4.84 and B^^Q = 
5.34, as inferred from solar photometric data from Allen 1973, § 75. and 
photometric transformations from Kent 1985 and Fukugita et al. 1995. 
Errors arc la. 

'^Notc that 0* and should likely be reduced by about 20%, due to the 
systematic effect discussed at the end of § 3. 



-17- 




-23 -22 -21 ^ -20 -19 -18 

- 5 log h 

Fig. 1. — 2(7 error ellipses in M* and a for the r-band luminosity functions from the CNOCl 
and Las Campanas redshift surveys. 
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-22 -21 -20 -19 -18 

Mb,„ - 5 log h, \ - 5 log h 

Fig. 2. — 2(7 error ellipses in M* and a for the 5-band luminosity functions of the CNOCl 
samples. Also shown are M* and a values for samples from the Canada-France (CFRS) and 
Autofib redshift surveys. Results for red and blue galaxy subsamples are shown only for 
CNOCl and CFRS. 
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-22 -21 -20 -19 -18 -17 

M, - 5 log h 

Fig. 3. — The r-band luminosity function of CNOCl samples. The lines are the STY fits, 
while the points are SWML solutions with la errors. Also shown is the STY fit from the 
Las Campanas survey. 
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Fig. 4. — The S^^-band luminosity functions of CNOCl samples are compared to those of 
corresponding CFRS samples. The points show SWML solutions for the CNOCl all sample. 
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Fig. 5. — The B-band luminosity function of the CNOCl all sample is compared to that of 
two samples from the Autofib redshift survey and to that of the z — 0.2 — 0.5 all sample 
from the CFRS. 
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Fig. 6. — The 5^B-band luminosity density pl for CNOCl samples, compared to that from 
analogous CFRS samples. The dashed hues show the CNOCl results after including a 20% 
reduction in normahzation as discussed in the text. Note that the CNOCl and CFRS curves 
differ primarily in their absolute normalization, particularly for the blue subsample. The 
shapes of the trends are very similar, once we rescale to match the normalizations. Also 
shown are results for local surveys as estimated by Lilly ct al. (1996), and in the top panel, 
the result from the LCRS, converted to the 5^B-band using Bab — Rlcrs — 1- 



